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ABSTRACT 


Extensions of the investigations aimed at preventing aircraft destruc- 
tion through the effects of hydraulic ram were conducted. By analyzing 
spark shadowgraphs the existence of a family of transient shock waves 
was confirmed. These shock waves were created by a surface anomaly in 
a projectile and are not unlike those found in fragmentary warheads. In an 
attempt to measure wall deflections subsequent to projectile impact, a 
system based on optical fiber proximity probes was investigated. Mathe- 
matical models based on grid distortions were used as a means of calculat- 
ing shock wave pressures. The AFIT thesis correlating cavity size, energy 
and radius to underwater explosions was verified from an analysis of 16 mm 


films. 
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I, INTRODUCTION 


The term hydraulic ram has been historically associated with the 
various waves causing, usually, catastrophic failure of an aircraft fuel 
cell when penetrated by a projectile. 

When a projectile or fragment with high velocity impacts the surface 
of a liquid medium, a disturbance is produced which Spread throughout 
the medium. The disturbance is characterized by a region of increased 
pressure and density together with a small forward motion of the fluid. 
When the energy of the disturbance is large, the disturbance travels as a 
surface of discontinuity and is called a shock wave. (Ref. 1). Typically 
the shock wave phenomenon has been identified as the first phase of hy- 
draulic ram. 

As projectile penetration continues, a pressure field is created that 
tends to-displace the fluid from the projectile, transforming projectile 
' energy to kinetic energy of the fuel. This is the beginning of the cavity 
phase. 

The phenomenon associated with these two phases has been studied 
for as long as three aes i.e. 1947 through present. Recent investigations 
into fragment-induced shock waves and cavities have beeninitiated, the 
results of which are leading to a more clear understanding of hydraulic 


ram. 
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II, SHOCK PHASE EXPERIMENTS 


A. WAVE GENERATION IN SHOCK PHASE BY FRAGMENTARY PROJECTILES 

Previous experiments conducted at NPSinvolving the impacting of 
0.22 caliber projectiles into a water-filled tank produced shadowgraphs 
such as those shown as Figures 1 and 2. Figure 1 is the projectile in free 
flight , whereas Figure 2 shows the family of waves created afterimpact. 
The flow features picible in Figure 2 are considered to be: 

1. projectile wake, 

2. reflected bandwave, 

3. reflected bow wave, 

4. incident wake recompression shock, 

5. incident bandwave, 

6. transient entry shock, 

7. secondary transient shock wave (see discussion below). 
In Figure 2, “items 6 to 9 occur in the water internal to the tank and items 
1 to 5 occur in the air external to the tank. The dense, thick region in 
each shadowgraph is the test tank wall. Region 8 is the projectile's 
Cavity and item 9 is the grid distortion used to determine the medium's 
density. The spark shadowgraph method was used as it has the advantage 
of showing the projectile-induced waves at one observation. Of particular 
interest was the formation of the wave 7, in Figure 2, which was investi- 


gated as follows. 
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0.22 Caliber Projectile in Free Flight 
Bow Shock Wave 
Band Shock Wave 


Wake Recompression Shock 
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A randomly selected, factory loaded 0.22 calprojectile has the geom- 
etry shown in Figure 3. The volume of the projectile's nose was calculated 
using ene eeS nine tmuncarcd cones while the portion aft of the shoulder was 
calculated assuming a solid cylinder. Assuming a constant velocity ene 
ately after penetration into the tank, the rate of volume penetration was 
plotted against distance along the projectile's centerline, Figure 4. The 
discontinuity in Figure 4 appears at a longitudinal distance of 0.20 inches 
and corresponds to the distance where the projectile itself has a discontin- 
uity, i.e., the shoulder. This shoulder is the origin of bandwave shown 
in Figure 1. This suggests that although the surface irregularity is small 
(0.020 inches), it is still capable of creating secondary transient shock 
waves. It is concluded that the wave 7 in question on Figure 2 is indeed 
a secondary transient entry shock. Reference 2 describes the deformation 
of several previous test rounds. Considering their fragmentary shape it 


could be expected that fuel cells punctured by fragmentation warheads 


would be subjected to a myriad of such secondary shocks. 


B. PRESSURE BY GRID DISTORTION 

Figure 5 is a shadowgraph of a fine-mesh grid taken while the fluid 
was atrest. No grid distortions are apparent. The projectile's path is 
from the left to right through the center of the shadowgraph. After penetra- 
tion the grid in this area is distorted. This region was used in the subse- 
quent analysis. 

As the projectile penetrates the water it displaces fluid by compressing 


it to very high pressures. This disturbance is the transient shock 6 in 
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Figure 2. As penetration continues, the transient shock continues to receive 
energy until the projectile becomes fully wetted; this occurs at 23 micro- 
seconds in Figure 4. A rarefaction wave then forms that tends to dissipate 
the shock wave. (Ref. 3). The localized high shock pressures are known 
to be of short duration, usually tens of microseconds, and are completely 
dissipated by the time significant projectile penetration occurs. 

The spark shadow-graph could be set off with any time delay up to 
approximately 30 microseconds. Each shock wave in the shadowgraph re- 
veals a distinct presence-of-light band and anabsence-of-light band. The 
concentration of light around the projectile's transient entry shock is best 
explained by considering the pressure and density increase due to the shock 
wave. Similar optical distortion at the nose of the fragment is due to high 
density at the nose. See Figure 2. The pressure is known to reach values 
in excess of 4000 kg/em2 (56,890 psi), which is sufficient to produce a 
change in the medium's refractive index. The width of the dark band varies 
with the pressure in the wave. The fractional change in refractive index, 
€ , has been studied by Poindexter and Rosen and was determined to vary 
asieo— il 06x 107°p where p is pressure in kg/em2 (Ref. 1). This conforms 
to an isothermal change, whereas a change produced by the transient shock 
is more nearly adiabatic (Ref. 1). A value of €= 107°p has been justified by 
McMillen and was used throughout this analysis (Ref. 1). It was also as- 
sumed that in a homogeneous medium such as water, light rays have the 


form of straight lines. 
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Fig. 6 Shock Front Pressure Generated in Water by Projectile Impact. 
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Seven shock pressures were extracted from Figure 6. When the pro- 
jectile had penetrated three cm into the tank the pressure was noted as 
2700 kg/em-. Define 8, and @2 as the angles which theincident ray and 


refracted ray make with a normal to the shock wave surface (Fig. 7). 





Fig. 7 Light rays in homogeneous regions 
bounded by a shock wave. 


By cae ee is given by (1) 
sin O,= 


uae em (2) 


Using Snell's law the angle Q1 can be introduced 
4 
Np smO,= NSM, 
sin 67 = (f \ser 8, 
Ny 


From the previous paragraph the change in refractive index can be expressed 


as €= AN = h- Ty 


TN 
mt (2) 
: Y= be 
From equations (1) and (2) 


sin 81=( (aes . 
“ Oy= sul’ ( LR (3) 


and ZA6 =(6,~ 81)+(°2-8/') = Z(8,- ,) (4) 
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Note that the changes in angle entering and leaving the shock wave are 
additive. If the density is homogeneous within the region bounded by the 
shock wave, as is assumed here, then the high pressure region behaves 
as a spherical lens. Forg = 0.90 from equation (1), 2) is found to be 
5.740°. For a pressure of 2700 kg/cm2, 8, is calculated to be 5.588° by 
equation (3); consequently 2A8 = 0.303°. Even for high pressures the 
angular deflections are small. Calculated values for the remaining six 
pressures are presented in Table 1. These tabulated results appear reason- 
able. 

The grid distortion analysis assumes rays have spherical symmetry 
i.e. the refractive index depends only on the distance r from a fixed point 


0; see Figure 8. 





6) 
Fig. 8 Illustrating Bouguer's formula nd=constant 


For rays situated in a plane through the origin, nrsin¢ =c, where cis 
a constant. With d=rsin$, nd=c 
This is known as the formula of Bouguer. From geometry, if (r,6) are 
the polar coordinates of a plane curve, the angle ) between the radius 


vector to a point P on the same curve and a tangent at P is given by 


surg = h(6) 


| r(8) 4 [ee i 
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Angular deflection of light ray by a spherical shock wave 


for various shock wave pressures 


Table 1 


Nondimensional radius 
Entry angle, radians 


RR 


Thetal 


Exit angle, radians 


Theta2 
Ansl 


Entry angle, degrees 


Exit angle, degrees 


2(ANS1-ANS2) 


Ans2 
Res 
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Angular deflection of light ray by a spherical shock wave 


for various shock wave pressures 


Table 1 


Nondimensional radius 
Entry angle, radians 
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Theta2 
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Angular deflection of light ray by a spherical shock wave 


for various shock wave pressures 


Table 1 


Nondimensional radius 
Entry angle, radians 
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It is possible to eliminate > by introducing the constant c, aS can be seen 


by the following manipulations: 
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The assumption of constant pressure behind the shock wave is not 
realistic. At this point, the assumption of uniform conditions behind the 
shock wave can be abandoned by introduction of equation (5). In equation 
(5), c is determined by the point where the ray intersects the shock wave. 
The index of refraction, n, which is relatedto pressure by E= 107°p, can 
be made a function of the radius. 

To adapt equation (5) to the specific case of the hydraulic ram shock 
phase itis Betsey to carefully define the various geometric parameters 
involved. Figure 9 illustrates the angles and radii necessary to apply 
equation (5). Note that the diagram is symmetric about r* at angle& . 

A functional form for n(r) can be assumed. For example, a parabolic 
relation would be 

n(r) = ar? + br+c 
In order to evaluate the constants a, b, and c, the integral of equation (5) 
must be-evaluated and iterated with geometric relations. This is quite, 
complicated . In lieu of a parabolic relation a linear formula ie used.'- This 
has the merit of a straightforward solution. One of the two parameters, i.e. 
A, can be determined leaving B to be found as the following manipulations 
show. 


Assume (Fig. 9) 


Wa np deal i Fy, Wet 
Vie 
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ah WSO ANG SS Vee (6) 
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From the law of Bouguer 

KV SUA $= constant 


WLYSUrG7 = ney* sug *= uhy* 
Nye sung = ‘a 


From analytic geometry and Snell's law 


SUR, = AL (8) 
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b= ae O'4e + 46 (10) 
= % +40-T - e'-2-26 + % (11) 
K= 18,4 LAB : (12) 
b> DT - 0'-g- ab +22, +246 | 

“BE q- QO +e 4+ 46 (13) 


ri , U 
From Figure 9 with cle NG) Ayr known 


am Aa a © visu B, 


Boe Mi c"sua( bo) = r/cun (B,- #) (14) 
Yaz 0 sua- v'sta(O-a) = 0 [suag! - sua(p'-x)] (15) 
fr +$,= + aQ 
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Using a pressure behind the shock wave of 400 kg/cm2, ng was 
found to be 1.335. See Figure 6 for data used in these calculations: 
r' is 7.62 cm. As shown in the Appendix the value for B is -0.04947. 
Using the equation for the linear pressure distribution it is calculated 
that p is 3646 kg/cm2 at & 0.1267. The pressure increases toward 
the projectile nose. This result is consistent with Figure 29 of Holm's 
thesis (Ref. 2). 

Results from fluid mechanics indicate that linear relations for p(r) 
and n(r) are not appropriate. Consequently the calculated pressure at 
4 = 0.1267 is probably not correct. This analysis gives the correct 


trends. The difference between experimental and calculated values is 


considered to result from the assumed distribution for n. 


C. WALL RESPONSE BY FIBER OPTICS PROXIMITY PROBES 

The structural response of the entrance wall during the shock phase 
has been investigated and is presented in Refs. 4,5, and 6. Since the shock 
wave is a SD sas disturbance propagating through the fluid medium, it is 
also by definition a stress wave. Types of stress waves include compression, 
tension and shear waves. The compression and tension waves involve posi- 
tive and negative pressures, respectively, where the media particle motion 
is in the direction of wave propagation. The shear wave involves particle 
motion that is normal to the wave direction. Compression and tension waves 
are also known as dilatation waves while shear waves are sometimes called 
distortion waves. Tension waves that reduce compressed areas are also 


known as rarefaction waves (Ref. 3). Loads transmitted through the tank's 
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structural members are initially transmitted as stress waves of both dilata- 
tion and shear form. The outward motion of the entrance wall somewhat 
relieves the pressure disturbances (stress waves) acting on the wall. To 
measure the outward displacement of the wall a system based on optical 
glass fiber bundles was constructed. The optic fiber bundles were arranged 
coaxially with the transmitting bundle on theinside. This light was then 
directed onto the aluminum entrance wall. The reflected light was collected 
and transmitted by the output fiber bundle to an RCA 931A photomultiplier 
tube and a Pacific Photometric Model 11 laboratory photometer. The fiber 
bundles were three feet in length and were comprised of three-mil flint 
glass fibers with an index of refraction of 1.62 and a cladding glass of in- 
dex 1.52. The corresponding numerical aperture was 0.56 with an accept- 
ance cone half angle of approximately 34 degrees. The glass bundles were 
secured in a mounting bracket and placed over the entrance wall on a radial 


line emanating from the projectile impact point. The mounting bracket al- 


“~ «Jowed-for:positioning the bundles up .to.a maximum of. nine inches from. the....0.., o 


impact point. To increase the intensity of reflected light the entrance wall 
was hand polished with #500 Wetordry Tri-M-ite paper, then with Crystalbay 
Crocus cloth. The signal output of the photomultiplier was transmitted to a 
dual beam oscilloscope fitted with a Polaroid camera. 

Initial calibration of the system was accomplished by mounting the 
fiber bundles in a platform micrometer. The light was reflected from a 
circular aluminum disk of the type used for the entrance wall. This disk 


was Similarly polished. A calibration curve of volts versus the distance 


separating the disk and glass bundles was obtained, Figure 10. 
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. Subsequent investigation of this system, modified to include a 
"chopper wheel", revealed an unacceptably low rise time for the photo- 
multiplier tube. A Pacific Photometric Model 62 multiplier tube (RCA model 
6199) was substituted; a Spectrum Physics model 132 HeNe laser replaced 
the light bulb. The inherent noisein the lattersystem obliterated the os- 
cilloscope trace, rendering this optical fiber concept useless without 


further development to improve signal-to-noise ratio. 
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Fig. 10 Calibration Curve for Optical Fiber Proximity Probe 
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III, CORRELATION OF OSCILLATION PERIOD, 


CAVITY SIZE, AND ENERGY DISSIPATION DURING CAVITY PHASE 


A. TECHNICAL APPROACH TO CORRELATION 

In displacing water from its path, the projectile imparts a radial 
velocity to the water. The momentum associated with this velocity field 
is sufficiently high to cause the flow to separate from the projectile and 
to continue to expand causing a cavity. Fluid pressures caused by the 
tank boundaries eventually halt the cavity's expansion. At this point 
most of the kinetic energy initially associated with the velocity field has 
been converted to stored energy, gravitational potential energy and de- 
formation of the tank's wall. These forces cause the cavity to start to 
collapse radially and, as the cavity shrinks, the fluid velocity increases 
resulting in the conversion of stored energy into kinetic energy and work 
in compressing gas trapped in the cavity. Since projectile impact was near 
the center of the tank the cavity would first neck down until anearly spheri- 
cal shape resulted. The sphere would then tend to collapse to a point until 
the pressures in the cavity halted this inward motion. Air which entered 
the cavity through the entrance wall hole is trapped and compressed to a 
high pressure at the instant of minimum cavity volume. A high pressure 
shock wave then propagates radially outward fromthe collapse point. After 
the pressure wave generated by the collapse subsides, the compressed air 
causes the cavity to re-expand. This cyclerepeats several times until all 
the stored energy is dissipated (Refs. 6 and 7). Various authors, Lundstrom, 
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Stepka and Cardea, et al., havelikened the cavity phase to underwater 
explosions. Using the method presented by Cardea and Cole, an analysis 
was made of a high speed film showing 0.30 cal and 14.5 mm rounds im- 
pacting a 6061T6 aluminum entrance wall. The film was obtained from 
Naval Weapons Center. These values were then compared with the analyti- 
cal solutions obtained by Holm's research at NPS. Two significant varia- 


tions were noted. Reference 7, equation (7), presents the maximum cavity 
B\i 
£o 


accurate. This difference accounts for an 8% variation in tabulated results 


radius as Oy SOT | A constant of -544 is considered to be more 
shown as Table 2. Secondly, this paper considers the effects of projectile 
impact through a pre-punched aluminum entrance wall whereas the 0. 30 

cal and 14.5 mm rounds impacted a solid tank wall. For the latter the 
shock wave is generated by the impact of the projectile on the wall vice 
the fluid medium. The generation and propagation will be identical ex- 
cept for wall density and sonic velocity effects. This problem is addressed 
- in Reference 3. The analysis of the 14.5 mm round must be considered as 
approximate as the maximum cavity radius to be meaSured extended beyond 


the frame size of 16 mm film. 


B. DATA OBTAINED FROM MOVIE FILM 

For each sequence of the 16 mm film obtained from the Naval Weapons 
Center, China Lake, the maximum cavity radius and time for one complete 
oscillation were measured. Using underwater explosion theory (Ref. 8) 


the energy required to expand the cavity was obtained first as a function 
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of cavity radius and second as a function of the cavity period. A theoret- 
ical value for the cavity radius was also calculated. These results were 


then compared with those obtained by Cardea (Ref. 7). 


C. RESULTS 

Table 2 lists the results obtained from the above measurements. The 
conclusion equating the cavity phase of hydraulicram to underwater ex- 
plosions appears valid. 

For a comparison of the time required for one complete cavity oscilla- 


tion, Reference 2 presents 


2 Poo \t 
is ¥\* 


where Oris the measured cavity radius. FromReference 8 this time is 


given as A n 
+= 4 Pot (37 Ve Gm 


Ss 
P, /e 
where 


Po= 2210 psf 
f5= 1.04 slugs 


Evaluating the constants yields 


nes Salhi Fe)" 
£= 1940 m [4e\" 


(~) 
respectively. It can be seen that the only difference in characteristic 


times is in the constants. Based on results of analysis of the film, 


equation is considered more correct. 
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0.30 cal 


Shot number millisec's 


nb whd 


14.5 mm 
Shot number 


aOunhwWDN HE 


t 


24.8 
22.0 
28.1 
28.8 
8) ol 


46.9 
48.3 
47.4 
44,4 
45.9 
43.9 


am 


measured 


Ft 


485 
2990 
632 
- 656 
. 632 


UG NAS 
1.094 
1.007 
1.234 
1.266 
1.234 


am 
calc 
Ft 


~455 
404 
016 
5 DAY) 
2934 


. 861 
. 887 
. 870 
~815 
- 836 
- 806 


Ec=f(Am) 
Ft-lbs 


874.78 
610.68 
1272.5 
1370.0 
1413.3 


5916.5 
6462.3 
GLOW sy, 
SOLS) 39 
9402.3 
4852.2 


Ec=£(T) 
Ft-lbs 


874.81 
610.70 
1272.6 
1370.1 
1413.3 


9916.7 
6462.5 


6107.7, 


9020.1 
9402.5 
4852.4 


Table 2 Cavity Period, Maximum Cavity Radius, Cavity Energy 


from Movie Film 


0.50 cal 


Shot Number millisec's 


on fw hd 


Table 3 Cavity Period, Maximum Cavity Radius, 


t 


24.8 
26.2 
A968) 
26.9 


am - 


observed 


eo 
anon on 


from Cardea (Ref. 7) 


am 
calc 


494 
aie 
-936 
-936 
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Ec 
Ft-lbs 


876.0 
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1117.0 
0 


Cavity Energy 
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VI. CONCLUSIONS AND PROPOSED RESEARCH 


A. CONCLUSIONS 

Secondary shocks will occur in fuel cells penetrated by fragmenta- 
- tion warheads. 

The mathematical model proposed for analyzing the optical grid dis- 
tortion yields correct trends. Further development of an assumed distri- 
bution for the index of refraction is warranted. 

Comparison of the cavity phase of hydraulic ram to underwater ex- 
plosions provides valid relationships between cavityradius, energy and 


period of oscillation. 


B. PROPOSED RESEARCH 

HIGH SPEED FRAMING CAMERA 

A high speed framing camera capable of 200,000 frames per second 
has been obtained. When incorporated into the esas test set-up, a 
complete time history of the shock front and projectile penetration will 
be available. 

STRAIN GAGES 

The existing entrance wall has been redesigned and fitted with several 
strain gages. Proposed research with this entrance wall should yield in- 
formation on wall deflections and wall bendingmoments caused by pro- 


jectile penetration. 
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APPENDIX A 


SAMPLE CALCULATIONS FOR THE GRID DISTORTION ANALYSIS 


line used for data 


Pen cy nA line through 
center 





Fig. 12 Enlarged Grid Distortion 
Steps for obtaining a solution: 
1. evaluate equation (17) to find Q, 
2. evaluate equations (8) and (9) to find Dey and Ab 
3. evaluate equation (12) to get & 
4. evaluate equations (10) and (13) for By and 9" 
5. evaluate equation (14) to get rt 
6. evaluate equation (15) to get Ur 


7. evaluate equation (16) for $, 


" 
L 


9. draw a curve of B vs r* using equation (18) 


8. evaluate equation (11) for 


10. using number pairs (B,r*) from step 9, evaluate equation (19) 


for the right hand side. Iterate until correct (B,r*) are found. 


36 











4 MRQMATIA 
Ss ~ { - a 
Ses NOImoree Ge a NOT sMOM AIOE, 
Pe | 
Blvd beau set! 


, 


fotef2i! bh osomlatd oft 
noaisoloe 2 malt 
i? aot) wi (l) oéHevoa saa 
ja Bas fy) boi 41 fb) bere (5) eno hfeuos Sai 

:. 
a = i} nov eyns sAsiae® 


" . | 
zi] brs 4 sak (24) ben (Ot) — " ojouwlave 


: 
iy dep oor (hb nebiedpe abeiiows 
4 90 on 29) relive ote 
a | 
~o iat fAljno po nape 
ad 
: ”“ 


vom GPa ioe 
7 (84) ; es oe i» week 
(C4 mownupe sAhlave ft! oy quit [4,8) y tone one 


sOnvo) S90 (99,2) Soares funy arial ole Virol tel ped tom 





From Figs. 11 and 12 
d = 0.1 inch 
L7 = 19.8 inches 
L, = 9 inches 
Yo = .35 inches 
[Pale GE 
P = 400 kg/cm? 
* A&A = (0.1) (180) = 0.299 = .005 radians 


(19. 8) 


From equation (17) 
0,35 = 35ua(6- ob) - (oe CDS (8,-~)| Tomo 
Assume (G-%) is small 
vs 0.352% 3(@,-A) - (9-3)% 
3&, = .35+ 9X 
- or Q, = 0.35 + 3% = .1318 radians 
3 


try 6,=1 55° 


§-a= 1,59°- 0.19= ues 
then 


0.35= 3 sur 11E-(9 - 30s 1,26} (+o. 0.29) 


from iteration 
ay 1.51° 
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Using p = 400 kg/em at r' = 3 inches 


é€= 10 p= 0.004 
W4= HK, ( [+4 .004) = 1.33 (1.004) 


Ro> (.3353 

aw Wye nse) 2. 0. 99CG6 
L 1.3365 

then 


Sun b= Eee <r) (.9960) 


Samby= (BIL 

6, (= ketoee 

Mb =) ib \= 1.57 - 1.5133" 
Eb= 0, 05\q 


from equation (12) 

= ote c 

z= .29- .051]7 = 0.0933" 

L 

from equation (10) 
e'= COT SMS eres, 90 t 0988 + OStT-1.5] 
6'= $2.515— 

from equation (13) 

6" = 90+2 +46-8,= 92.515" 

from equation (14) 
Wi= SUC -2)- Y= 35 (1.57- 119) — .035 
Y= 0,03" 
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from equation (15) 
yas B [sue TST Sum (1.9- 19)| 
Y= 0.015" 


from equation (16) 


Gi= 190 + .0>-1,51 


Ol = [th a 


from equation (11) 


a 
fe hae =e] Osh 
. = 1S 
An approximation for B is 
Kk 
Ce et 
mie 2 Uji = O14 2-03 2 (7L7 
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B= (1.3363). 1508 1261) 
(.(261)Cl- 1267) 


B=— 0.04947° 


from the assumed distribution for n 


w= Wy + ie 
r! 


R= 1.3353 + 0.04947 (I-.1t6T) when 
N= ls Ono> 
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